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Abstract 

We present a measurement of the K-tt vector current form-factor parameters for 
the decay K L -> ir^e^v. We use 328 pb" 1 of data collected in 2001 and 2002, 
corresponding to ~ 2 million K e s events. Measurements of semileptonic form factors 
provide information about the dynamics of the strong interaction and are necessary 
for evaluation of the phase-space integral If^ needed to measure the CKM matrix 
element \V US \ for K L -> n^e^v decays. Our result is A + = (28.6 ± 0.5 ± 0.4)10~ 3 
for a linear fit, and X' + = (25.5 ± 1.5 ± 1.0)10" 3 , A'| = (1.4 ± 0.7 ± 0.4)10~ 3 for a 
quadratic fit. 

key words: ke3 form factor 
PACS: 13.20.Eb 
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1 Introduction 



Semileptonic kaon decays, Kl — > n ± £ :f h', (Fig. 1) offer possibly the cleanest 
way to obtain an accurate value of the Cabibbo angle, or better, V us . Since 
K — » 7r is a 0" — > 0~ transition, only the vector part of the weak current 
has a nonvanishing contribution. The transition is therefore protected by the 
Ademollo-Gatto theorem against SU(3) breaking corrections to lowest order. 
At present, the largest uncertainty in calculating V us from the decay rate, 
is due to the difficulties in computing the matrix element (n\j^\K). In the 




Fig. 1. Amplitude for Kl — ► n i^u. The gray region indicates the K — ► ttW vertex 
structure. 

electron mode Kl — > 7r ± e =F z/ only one form factor is involved. In the following 
we will use the notation shown in Fig. 1, in which P, p, k and k! are the 
kaon, pion, electron and neutrino momenta, respectively; m is the mass of the 
charged pion and M that of the neutral kaon. Terms in (P — p) M that acquire 
factor of m e are neglected. Therefore: 

Wl\K) = f + (t){P + p), 

We replace the form factor above with /+(0) f+{t), where t = (P — p) 2 = 
(k + k') 2 = M 2 + m 2 — 2ME n is the only L-invariant variable and /+(0) = 1. 
The form factor is dominated by the vector K-ir resonances, the closest being 
the if* (892). Note that for t > 0, f + (t) > 1. The presence of the form factor 
increases the value of the phase-space integral and the decay rate. The natural 
form for f(t) is 

It is also customary to expand the form factor as 

hit) = i + K— 2 + \K ( 2 ) 

m 2 2 \mv 

In the following we retain linear and quadratic terms. Note that the expansion 
of the pole form above gives A'+ = (m/M v ) 2 and A'j = 2 A' + . From P + p = 
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2P — [k + k'), and neglecting the k + k 1 term which is also proportional to m e , 
the amplitude is: 

Tt = (neu\H w \K) = 2 WtM ^(^(i - ^)u{k>) f + {t) . (3) 



Squaring, summing over spins, and integrating over all variables but the pion 
energy, we obtain the pion spectrum 

g{z)<x (, 2 -4a) f (l + ^(£-z) + ^(£-z) 2 V 

where a = m 2 /M 2 , £ = 1 + a and z = 2 E n /M is the normalized pion energy. 
The spectrum can also be written in terms of t = m 2 (£ — z)/a as 



h(t) oc \ i/2 f + (t) 2 

\ = t 2 -2t(M 2 + 2m 2 ) + (M 2 -m 2 ) 2 



(4) 



2 The KLOE detector 



The KLOE detector consists of a large, cylindrical drift chamber (DC), sur- 
rounded by a lead/scintillating-fiber electromagnetic calorimeter (EMC). A 
superconducting coil around the calorimeter provides a 0.52 T field. The 
drift chamber [1] is 4 m in diameter and 3.3 m long. The momentum res- 
olution is cr p± /p± ~ 0.4%. Two-track vertices are reconstructed with a spa- 
tial resolution of ~ 3 mm. The calorimeter [2] is divided into a barrel and 
two endcaps. It covers 98% of the solid angle. Cells close in time and space 
are grouped into calorimeter clusters. The energy and time resolutions are 
a E /E = 5.7%/^jE (GeV) and a t = 57 ps/y^ (GeV) © 100 ps, respectively. 
The KLOE trigger [3] uses calorimeter and chamber information. For this 
analysis, only the calorimeter signals are used. Two energy deposition above 
threshold (E > 50 MeV for the barrel and E > 150 MeV for the endcaps) are 
required. Recognition and rejection of cosmic-ray events is also performed at 
the trigger level. Events with two energy deposition above a 30 MeV threshold 
in the outermost calorimeter plane are rejected. 



3 Analysis 



The 328 pb" 1 of 2001-2002 data used in this analysis [4], is divided into 14 
periods of about 25 pb -1 /period. For each data period we have a corresponding 
sample of Monte Carlo events with approximately the same statistics. 
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Candidate K L events are tagged by the presence of a K s — > ix + ix~ decay. The 
Kl tagging algorithm is fully described in [5] and [6]. The Kl momentum, 
Pk l , is obtained from the kinematics of the <ft — > K$Kl decay, using the 
Kg direction reconstructed from the measured momenta of the decay tracks 
and the known value of p^. The resolution is dominated by the beam-energy 
spread, and amounts to about 0.8 MeV/c. The position of the cj) production 
point, x,£, is determined as the point of closest approach of the Ks, propagated 
backward from the Ks vertex, to the beam line. The Kl line of flight {tagging 
line) is then constructed from the Kl momentum, px L = P</> — Px s , and the 
position of the production vertex, x^. 

The efficiency of the tagging procedure depends slightly on the evolution of 
the Kl, mainly because the trigger efficiency depends on the Kl behavior. 
To identify events in which the Ks by itself satisfies the calorimeter trigger, 
we require the presence of two clusters from the Ks — > 7r + 7r~ decay associated 
with fired trigger sectors (autotrigger). The value of the tagging efficiency 
obtained from Monte Carlo is about 40 % and is independent of t to within 
0.4%. 

All tracks in the chamber, after removal of those from the Ks decay and their 
descendants, are extrapolated to their points of closest approach to the tagging 
line. For each track candidate, we evaluate the point of closest approach to the 
tagging line, x c , and the distance of closest approach, d c . The momentum p c 
of the track at x c and the extrapolation length l c are also computed. Tracks 
satisfying d c < ar xy + b, with a = 0.03 and 6 = 3 cm, and —20 cm < l c < 25 
cm are accepted as Kl decay products, where r xy is the distance of the vertex 
from the origin in the transverse plane. For each sign of charge we consider the 
track with the smallest value of d c to be associated to the Kl decay. Starting 
from these track candidates a vertex is reconstructed. The combined tracking 
and vertexing efficiency for K e ^ decays is about 54%. It is determined from 
data as described in Ref. 5. An event is retained if the vertex is in the fiducial 
volume 35 < r xy < 150 cm and \z\ < 120 cm. 

To remove background from Kl — > 7r + 7r _ 7r° and Kl — * decays with min- 

imal efficiency loss, we apply loose kinematic cuts: assuming the two tracks 
to have the pion mass, we require E^ iss — "p 2 miss — M% < —5000 MeV 2 and 

\J Emiss + Pmiss > 10 MeV, where E miss and p m iss are the missing energy and 
momentum, respectively. A large amount of background from Kl — > ~k\w de- 
cays is rejected using A^, the lesser value of \E miss — p m i SS \ calculated in the 
two hypotheses, nfi or fin. We retain events only if this variable is greater 
than 10 MeV. After the kinematic cuts described above, the efficiency for the 
signal is about 96%. 

These kinematic criteria do not provide enaugh suppression of the back- 
ground from Kl — > vr ± e =F z/ decays with incorrect charge assignment and from 
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Kl — > TTfiu decays. We make use of time-of-flight (TOF) information to further 
reduce the contamination. 



For the purpose of track-to-cluster association, we define two quantities related 
to the distance between the track, extrapolated to the entry point of the 
calorimeter, and the closest cluster: g?tc ; the distance from the extrapolated 
entry point on the calorimeter to the cluster centroid and drc±, the component 
of this distance in the plane transverse to the momentum of the track at the 
entry position. We only consider clusters with dq<c± < 30 cm and E >50 MeV. 

We evaluate the cluster efficiency using the Monte Carlo, and correct it with 
the ratio of data and Monte Carlo efficiencies obtained from control samples. A 
sample of K e % events with a purity of 99.5% is selected by means of kinematics 
and independent calorimeter information. Figure 2 shows the corrections as 
a function of t obtained for a single run period. It is worth emphasizing that 

1.02 



1.00 - 



0.98 - 



0.96 - 



0.94 



0.92 



0.90 
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Fig. 2. Cluster efficiency correction for pions (ir in red, ir + in black) on the end 
cap from control sample. The plot refers to a single period of data taking 

if this correction were not taken into account, the effect on A' + would be 
large (about 20%) and would produce different results for each charge (about 
15%). For this reason, the analysis is performed separately for each charge. 
The comparison of the two results provides a first check of the validity of the 
corrections. 

For each Kl decay track with an associated cluster, we define the variable: 
Atj = t c i — ti, (i = 7r, e) in which t c i is the cluster time and U is the expected 
time of flight, evaluated according to a well-defined mass hypothesis. The 
evaluating of ti includes the propagation time from the entry point to the 
cluster centroid [7]. We determine the e + e~ collision time, t , using the clusters 
from the K$- 
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An effective way to select the correct mass assignment, Tie or eir, is obtained 
by choosing the lesser of \At n + — At e -| and \At n - — At e +\. After the mass 
assignment has been made, we consider the variables Atn + At e and At n — At e . 
These are shown in Fig. 3 for signal and background Monte Carlo events. We 
select the signal by using a 2a cut, where the resolution a ~ 0.5 ns. After the 



At e + At, (ns) 




Fig. 3. Monte Carlo distribution of Aie + At^ versus At e — At n . Signal (gray scale), 
background from opposite sign Kl — > ir^e^v (black dot), background from same 
sign Kl — > npiv (empty box). 

TOF cut we have a contamination of ~0.7%, almost entirely due to Kl — > tiliv 
decays. 

We take the TOF efficiency from the Monte Carlo after correcting the time 
response of the calorimeter using data control samples [7]. The quality of this 
correction can be checked by comparing the data and Monte Carlo distribu- 
tions for At e and At n shown in Fig. 4, obtained from the same control sample 
used for the cluster efficiency. 

We measure the form-factor slope parameters by fitting the distribution of 
the selected events in t/m 2 . We modify the kinematic range of t/m 2 , varying 
from (m e /m) 2 ~ up to (M — m) 2 /m 2 ~ 6.8, to [-0.5,6], to take into 

account the smearing effect at t ~ and the low statistics at high values of 
t. After subtracting the residual background as estimated from Monte Carlo, 
we perform the fit using the following formula: 

dN 20 

— (i) = N J^A(i,j) x p(j,X' + ,K) x e tot (j) x F FSR (j) (5) 
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(a) (b) 
Fig. 4. Atj for electron (a) and pion (b) for data (black) and Monte Carlo (red) 

where p(j, A' + ,A'|) is the three-body differential decay width as defined in 
Eq. (4), and A(i,j) is the probability that an event with true value of t/m 2 
in the j th bin has a reconstructed value in the i th bin. The chosen bin size is 
0.5, which corresponds to about 1.6 a t , where a t is the resolution in t/m 2 . 

The total efficiency, e tot (t), takes into account the acceptance and the effi- 
ciency of the analysis cuts. -Ffsr is the correction for final-state radiation. It is 
evaluated using the KLOE Monte Carlo simulation, GEANFI [8] , where FSR 
processes are simulated according the procedures described in Ref. 9. FSR af- 
fects the t-distribution mainly for high-energy pions, i. e. for low t, where the 
correction is 3-5%. The slopes A' + and A" are free parameters in the fit while 
the No constant is the total number of signal events. 



4 Systematic uncertainties 



The systematic errors due to the evaluation of corrections, data-Monte Carlo 
inconsistencies, result stability, momentum miscalibration, and background 
contamination are summarized in Table 1. 

We evaluate the systematic uncertainty of the tagging efficiency by repeating 
the measurement using a tagging algorithm without the requirement of the 
autotrigger. We observe a change of 0.14xl0 -3 for A + in the case of the linear 
fit, and changes of 0.18xl0 -3 and 0.02 xl0~ 3 for A' + and A" respectively for 
the quadratic fit. 

We evaluate the systematic uncertainties on the tracking efficiency corrections 
by checking stability of the result when the track selection criteria are modi- 
fied. We establish the validity of the method by comparing the efficiencies from 



6 





Linear fit 


Quadratic fit 


Source 


5\+ x 10 3 


5\' + x 10 3 


5\'[ x 10 3 


Tagging 


0.14 


0.18 


0.02 


Tracking and vertexing 


0.16 


0.22 


0.18 


Clustering 


0.07 


0.24 


0.13 


Time-of-flight 


0.29 


0.87 


0.27 


Background 


0.08 


0.16 


0.03 


Momentum-scale 


0.06 


0.05 


0.05 


Momentum-resolution 


0.17 


0.22 


0.19 


Total systematic 


0.42 


0.98 


0.39 



Table 1 

Summary of systematic uncertainties on A+, A' + and A". 

data and Monte Carlo control samples, and from the Monte Carlo truth [7]. 
The uncertainty on the tracking efficiency correction is dominated by sample 
statistics and by the variation of the results observed using different criteria 
to identify tracks from Kl decays. The correction is run-period dependent; 
its statistical error is taken into account in the fit. We study the effect of dif- 
ferences in the resolution with which the variable d c is reconstructed in data 
and in Monte Carlo events, and the possible bias introduced in the selection 
of the control sample, by varying the values of the cuts made on this variable 
when associating tracks to Kl vertices. For each variation, corresponding to 
a maximal change of the tracking efficiency of about ±15%, we evaluate the 
complete tracking-efficiency correction and measure the slope parameters. We 
observe a change of 0.16xl0~ 3 for A + in the case of the linear fit, and changes 
of 0.22xl0 -3 and 0.18xl0 -3 for A' + and A" , respectively, for the quadratic 
fit. We find a smaller uncertainty by comparing the efficiencies from data and 
Monte Carlo control samples, and Monte Carlo truth. However, we conserva- 
tively assume the systematic uncertainty to be given by the changes in the 
result observed by varying the cut on d c . 

We evaluate the systematic uncertainties on the clustering efficiency correc- 
tions by checking stability of the result when the track-to-cluster association 
criteria are modified. In this case as well, the uncertainties on the clustering ef- 
ficiency corrections are dominated by sample statistics and by the variation of 
the results observed using different criteria for the track-to-cluster association. 
The correction is run-period dependent; we take into account its statistical er- 
ror in the fit. The most effective variable in the definition of track-to-cluster 
association is the transverse distance, g?tc±- We vary the cut on dr C ± in a wide 
range from 7 cm to 30 cm, corresponding to a change in efficiency of about 
17%. For each configuration, we obtain the complete track extrapolation and 
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clustering efficiency correction and use it to evaluate the slopes. We observe 
a variation of 0.07xl0 -3 for A + in the case of the linear fit, and variations 
of 0.24xl0 -3 and 0.13xl0 -3 for A' + and A" , respectively, in the case of the 
quadratic fit. We find a comparable uncertainty for A+ and for A" by com- 
paring the efficiencies from data and Monte Carlo control samples, and the 
Monte Carlo truth. 

We study the uncertainty on the Monte Carlo efficiency of the TOF selection 
procedure by measuring it using a pure K e3 control sample, and using the ratio 
of data and Monte Carlo efficiencies estimated in this way as a correction. 1 
The control sample is selected using tighter kinematic cuts and the calorimeter 
particle identification described in Ref. 5. The contamination of the control 
sample amounts to 0.4%. When applying the correction, we find a change 
in the result of 0.29 xl0~ 3 for A + in the case of the linear fit, and changes 
of 0.87xl0~ 3 and 0.27xl0~ 3 for A' + and A" , respectively, in the case of the 
quadratic fit. These variations are well within the statistical uncertainties. 

We study the uncertainties on the background evaluation by repeating the 
measurement on a sample with reduced background contamination. This is 
achived by identifying the electron using the energy deposition in the calorime- 
ter combined with a neural network (NN). In Fig. 5 we show the distribution 
of the NN output for the sample used in the analysis. Requiring the value of 
the NN output to be greater than 0.4, we reduce the background contamina- 
tion by about a factor of three. The differences in the result obtained with 
this cut are 0.08xl0~ 3 for A + in the case of the linear fit, and 0.16xl0 -3 and 
0.03 xl0~ 3 for A' + and A" , respectively, in the case of the quadratic fit. 

The effect of the momentum scale uncertainty and the momentum resolution 
have also been considered. We find the following relations by changing the 
momentum scale: 

Jp 5\>l bp 

a; - P ' a'; " P 

We conservatively assume a momentum scale uncertainty of 0.1%, which is 
much greater than the value obtained from a dedicated analysis [10]. This 
translates into a change of 0.06xl0 -3 for A + in the case of the linear fit, and 
changes of 0.05 xl0~ 3 and 0.05 xl0~ 3 for A' + and A" , respectively, in the case 
of the quadratic fit. 

We investigate the effect of the momentum resolution by changing the value of 
the resolution on t/m 2 . A variation of 3% worsens the fit quality, giving a x 2 



The TOF corrections cannot be used directly in the analysis because of the cor- 
relation between the energy response in the calorimeter and the TOF. 



8 



10 4 



10 3 



10 2 



10 



-0.5 0.5 1 1.5 

NN output 

Fig. 5. NN output distribution for data (dots), Monte Carlo signal (black solid line) 
and background (red solid line). 

probability variation of one standard deviation. The corresponding absolute 
changes are 0.17xl0~ 3 for A + in the case of the linear fit, and 0.22xl0 -3 and 
0.19xl0 -3 for X' + and A" in the case of quadratic fit. Varying the resolution 
of t/m 2 by a larger amount (6%) gives an unacceptable x 2 probability, about 
10 -9 , while nearly the same variations for the fit parameters are observed. In 
principle, if the distribution has a linear behavior, the slope is insensitive to 
any smearing due to the resolution. The only effect is due to the depletion 
of the bins at the boundary of the distribution, which worsens the \ 2 °f the 
fit. We have verified that the sensitivity to the momentum resolution is much 
smaller for a reduced fit range. 



5 Results 

About 2 million K e3 events were selected. The results of the linear fit obtained 
from all run periods are given in Table 2. The fit is performed separately for 
Kl — > 7T + e~V and — > n~e + v events to check the reliability of the evalua- 
tion of the cluster efficiency The results are consistent only if the respective 
efficiency corrections for each pion charge are applied. Then, combining the 
two charge results and including the systematic uncertainties listed in Table 1 
we obtain: 

A+ = (28.6 ± 0.5 stat . ± 0.4 syst .) x 1(T 3 
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The results obtained for the quadratic fit are given in Table 3. A correlation 
of ~ —0.95 between the A' + and A" parameters is obtained, as expected from 
the form of the parametrization in Eq. (2). A very slight preference for a 
small quadratic term is observed as indicated by the small improvement in 
the fit probability going from the linear, P(x 2 ) = 89%, to the quadratic fit 
P(x 2 ) = 92%. Including the systematic uncertainties listed in Table 1 we 
obtain: 

A' + = (25.5 ± 1.5 stat . ± 1.0 syst .) x 10~ 3 
A'| = (1.4±0.7 stat .±0.4 syst .) x 10- 3 

Figure 6 shows the t/m 2 distribution for the data and the fit result. The ratio 
data/fit is also shown. 



Linear fit 


A+ x 10 3 


X 2 /ndf 


Kl — > TT~e + v 
Kl — > Ti + e~V 
Combined 


28.7 ±0.7 

28.5 ±0.6 

28.6 ±0.5 


156/181 
174/181 
330/363 



Table 2 

Fit results in the hypothesis A" =0. Only statistical errors are shown. 



Quadratic fit 


A' + x 10 3 


A'| x 10 3 


X 2 /ndf 


Kl — ► 7r~ e + u 
Kl — > ir + e~V 
Combined 


24.6 ±2.1 

26.4 ± 2.1 

25.5 ± 1.5 


1.9 ± 1.0 
1.0 ± 1.0 
1.4 ±0.7 


152/180 
173/180 
325/362 



Table 3 

Fit results with A' + and A'_{_ as free parameters. Only statistical errors are shown. 

We also fit the data using the one-pole parametrization (see Eq. (1)). We 
obtain M v = (870 ± 6) MeV with x 2 /ndf = 326/364 and a probability of 
P(x 2 ) = 92.4%. Taking the systematic error into account, we obtain: 

M 1/ = (870±6 stat .±7 syst .)MeV 

This result indicates that, although the pole is dominated by the K* vector 
meson, contributions from other J p = 1~ resonant and non- resonant Kir 
scattering amplitudes are not negligible. 
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Fig. 6. Fit results: data (dots) are superimposed on the fit function (histogram). 
The data/Fit ratio is also shown. 



Conclusion 



We have obtained precise new values of the slopes used to describe the hadronic 
form factor in K e3 decay. The new KLOE result is consistent with the presence 
of a small quadratic term in agreement with the expectation of the one-pole 
expansion. The value of A + obtained with the linear fit is in good agreement 
with other existing measurements. The comparison with other existing mea- 
surements is shown in Fig. 7 in the case of the quadratic fit. Our result is in 
good agreement with ISTRA+ [11] and NA48 [12] and in marginal disagree- 
ment with KTeV [13]. 
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and Part. Astrophys.' of Deutsche Forschungsgemeinschaft, Contract No. GK 
742; by INTAS, contracts 96-624, 99-37; by TARI, contract HPRI-CT-1999- 
00088. 
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